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X-ray absorption by matter has long been described by the famous Beer-Lambert law. Here we 
show how this fundamental law needs to be modified for high-intensity coherent x-ray pulses, now 
available at x-ray free electron lasers, due to the onset of stimulated elastic forward scattering. We 
present an analytical expression for the modified polarization-dependent Beer-Lambert law for the 
case of resonant core-to-valence electronic transitions and incident transform limited x-ray pulses. 

Upon transmission through a solid, the absorption and dichroic contrasts are found to vanish with 
increasing x-ray intensity, with the stimulation threshold lowered by orders of magnitude through 
a super-radiative coherent effect. Our results have broad implications for the study of matter with 
x-ray lasers. 


Non-linear interactions of intense electromagnetic ra¬ 
diation with matter have long been utilized in the mi¬ 
crowave and optical regions to control nuclear and va¬ 
lence electronic transitions and have enabled break¬ 
throughs in many fields of science, such as medical imag¬ 
ing, telecommunication or the creation and manipulation 
of novel states of matter. The natural extension of these 
techniques into the x-ray region had to await the avail¬ 
ability of sufficiently bright x-ray sources in the form of 
x-ray free electron lasers. Over the last few years, several 
experiments performed with rather uncontrolled x-ray 
pulses of high intensity, produced through the self am¬ 
plification of spontaneous emission (SASE) process Q}, 
have revealed the presence of high intensity effects due 
to electronic stimulation Q or multiple ionization Q. 

Here we discuss how x-ray transmission through mat¬ 
ter can be modified in a controlled way by stimulated 
scattering effects induced by transform limited x-ray 
pulses now available through self-seeding Q. In contrast 
to stimulated inelastic scattering [2j,|9|,|l0[, which requires 
pulses with a broad bandwidth that covers the difference 
between excitation and de-excitation energies or multi¬ 
color pulses with separate “pump” and “dump” func¬ 
tions, we consider here the conceptually simpler case of 
elastic stimulation which exists within the energy band¬ 
width of the incident beam itself. In this case stimulated 
x-ray scattering modifies the fundamental Beer-Lambert 
law because of the direct link of x-ray absorption and 
resonant elastic scattering through the optical theorem. 

Of particular importance and interest are experiments 
that utilize resonant electronic core-to-valence transi¬ 
tions since they exhibit large cross sections, and for solids 
provide elemental and chemical bonding specificity, and 
through their polarization dependence enable the deter¬ 
mination of bond orientation [5[ and the dichroic separa¬ 
tion of charge and spin based phenomena @|. Resonant 
x-rays are widely utilized in x-ray absorption, x-ra y s cat- 
tering and coherent x-ray imaging experiments [5|-[8| . 

We derive the modified Beer-Lambert law by utilizing 
the time-dependent density matrix approach where the 
evolution of the resonant core-valence two-level system is 


governed by the optical Bloch equations [LI] • An analyti¬ 
cal solution is obtained for the case of incident transform 
limited x-ray pulses whose coherence time is much longer 
than the core hole lifetime. We apply our theory to the 
important case of 3d transition metal samples whose po¬ 
larization dependent transmission exhibits both a charge 
and spin response, the latter through the x-ray magnetic 
circular dichroism (XMCD) effect. We find that for the 
prominent Co L3 absorption resonance at 778 eV (wave¬ 
length of 1.6 nm), stimulated decays begin to rob inten¬ 
sity from the dominant spontaneous Auger channel at 
an incident intensity of about lmJ/cm 2 /fs (lTW/cm 2 ), 
with the onset lowered by a super-radiative coherent scat¬ 
tering enhancement in the forward direction. At higher 
intensities the sample becomes increasingly transparent 
with the spin-based XMCD contrast disappearing sooner 
than the charge-based absorption contrast. 

We follow the formalism of reference [fij and, denoting 
the x-ray polarization by the labels q = 0 for linear, q = + 
for right and q = — for left circular polarization, describe 
the polarization dependent x-ray response of a magnetic 
sample in terms of the atomic scattering length in the 
soft-x-ray approximation as f q (Q = 0) = r 0 Z+f ,q — i f" q , 
where ro is the Thomson scattering length and Z the 
atomic number. The spontaneously transmitted inten¬ 
sity through a sample of atomic number density p a and 
thickness d is given by the Beer-Lambert law 


Tl = Tl „-2A/"Vad 

^trans e 


(i) 


where I q rans and /g are the polarization dependent trans¬ 
mitted and incident intensities and oA = 2A f" q is the x- 
ray absorption cross section. Our x-ray scattering length 
formulation is related to the optical constants and the 
electric susceptibility through the complex refractive in¬ 
dex n q = 1 — 6 q + i/3 9 ~ 1 + r;(x /q + ix" q )> where 
S q = p a \ 2 (r 0 Z + f ,q )/2TT and p q = p a \ 2 f" q /2ir. The 
resonant polarization dependent x-ray absorption cross 
section a ahs = 2\f" q and the differential atomic elas¬ 
tic scattering cross section da q cat /dUl = ( f ,q ) 2 + ( f " q ) 2 
have a Lorentzian lineshape and are linked by the optical 
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theorem which may be written as, 

fiiq_ r ^7T r. , q . 2 . firq\21 _ (T/2) 2 _ /<j\ 

J Tl A LU j J J T 27r(^-£ 0 ) 2 +(r/2) 21 ' J 

Here £o is the resonant photon energy, T = + Ta is 

the total spontaneous decay width, which in the soft x- 
ray region is dominated by the Auger width r ~ Ta [13 ■ 
The polarization dependent radiative transition widths 
T| consist of a radial and angular part and can be calcu¬ 
lated by ab initio methods. We have derived their values 
for Fe, Co and Ni metal from experimental data, and 
they are listed in TableU]. 


TABLE I: Polarization dependent parameters for the L 3 res¬ 
onances of Fe, Co and Ni metals. Listed are the atomic num¬ 
ber densities p a , the resonance energies and wavelengths, and 
the polarization dependent (g = 0 , ±) peak experimental cross 
sections CTq, assuming propagation along the magnetization 
direction. F^ is the polarization dependent dipole transition 
width which includes the number of valence holes Ah, and F 
is the natural decay energy width [T3|. 
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Fe 

84.9 

707 

1.75 

8.8 

6.9 

5.0 

1.37 

1.08 

0.78 

0.36 

Co 

90.9 

778 

1.59 

7.9 

6.25 

4.65 

1.208 

0.96 

0.715 0.43 

Ni 

91.4 

853 

1.45 

5.1 

4.4 

3.7 

0.675 

0.575 

0.48 

0.48 


The polarization dependent Lorentzian x-ray absorp¬ 
tion cross sections cr^ bs = 2A f" q calculated with Eq.[2] 
and the parameters for Co in TableU] are shown as blue 
curves in Fig.Q](a). They were derived from fits of the 
experimental resonant cross sections by Voigt profiles as 
shown in Fig.Q](b), consisting of a convolution of the 
natural Lorentzian lineshapes in (a) with a Gaussian of 
1.4 eV FWHM to account for the band-structure broad¬ 
ened d valence states into which the 2p 3 / 2 core electrons 
are excited. 

For a sample of finite thickness d and atomic num¬ 
ber density p a , the transmitted intensity decays exponen¬ 
tially with the number of atoms in the beam N a /A = p a d 
according to Eq.Q] Since absorption and resonant scat¬ 
tering are related through Eq.O the Beer-Lambert ab¬ 
sorption law can also be derived by considering resonant 
elastic forward scattering. To do so one considers scat¬ 
tering by a thin atomic sheet so that the first Born ap¬ 
proximation is valid. For a sheet thickness A <C A the 
spontaneously forward scattered fields are coherent and 
the transmitted field is given by, 

ELns = El e ikA {1 - i A [r 0 Z + f q -i f" q ] Pa A} (3) 

Neglecting the non-resonant term tqZ , the intensity 
transmitted through the sample with N a atoms in the 
beam of cross sectional area A is, 

l^ran S | 2 =l^| 2 {l-2A^/" 9 +A 2 ^|[(/ /9 ) 2 +(/" 9 ) 2 ]J(4) 
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FIG. 1: (a) Polarization and photon energy dependent L 3 ab¬ 
sorption cross sections < 7 ^ bs = 2 A f" q for Co metal, calculated 
by use of Eq.Q] and the parameters T and in Table]!] (b) 
Comparison of the experimental dichroic cross sections (black 
lines) and the theoretical cross sections in (a) convoluted with 
a Gaussian of 1.4 eV FWHM. Note that the corresponding 
blue and red curves have the same areas but their peak val¬ 
ues differ by a factor of <Tj/ctq = 2.9. 


The first term is the incident intensity and the second 
term is the absorption loss (minus sign) in linear re¬ 
sponse. Within the Born approximation, the absorp¬ 
tion loss arises from the destructive interference of the 
incident field with the coherently forward scattered field. 
The third term is the forward scattered gain (plus sign) 
due to the coherent superposition of the fields scattered 
by the atoms in the sheet which scales with N 2 . In Eq.Q] 
we have neglected the weak intensity (8ir/3)N a [(f' q ) 2 + 
i.f" q ) 2 ]/A which is incoherently scattered. The coherent 
forward scattered intensity is larger by the enhancement 
factor, 


Qcoh — 0 N a 

»7r A 


( 5 ) 


where dfl C oh = A 2 /A represents the solid angle of co¬ 
herent forward scattering. The total field transmitted 
through a sample of arbitrary thickness d = NA is 
obtained by using the Darwin-Prins dynamical scatter¬ 
ing formalism to sum Eq.Q] over N thin sheets, which 
yields the exponential Beer-Lambert law Eq.Q] 13]. Re¬ 
markably, for forward scattering, the longitudinal co¬ 
herence length £ c ~ A 2 /AA does not enter [3]) and 
both Eq.Q] and Q c D h depend on the number of atoms 
N a in the beam cross sectional area A, and not on the 
sample thickness d. Q C oh is thus the same for a thin 
film with p a — 100atoms/nm 3 and typical thickness 
d ~ l/(cr a bs/9 a ) — 20 nm and a gas sample of the same 
atoms of density p a — 0.01 atoms/nm 3 and a much larger 
thickness of d ~ 200 pm. 

As the incident intensity is increased, the Kramers- 
Heisenberg perturbation theory (see 0) leads to un¬ 
physical results since it does not account for popula¬ 
tion changes in the excited state. This is overcome by 
the density matrix formalism which directly calculates 
the time-dependent ground and excited state populations 
[3 , which we shall denote pu(t) and P 22 (t) = 1 — pn(t), 
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respectively. The populations are obtained as the solu¬ 
tions of the optical Bloch equations. 

In the presence of stimulation, we can write the atomic 
scattering length as the sum of a spontaneous (subscript 
“0”) and stimulated non-linear (subscript “NL”) part ac¬ 
cording to, 

r = r 0 z+n +n L , r=n +n L (6) 

and Eq.Q] is replaced by, 



The spontaneous absorption cross section er a b s = 2A f" q 
with f' q = /q q in Eq.Q] becomes the coherence time and 
intensity dependent expression, 

<r abs = 2A [f" q 0 + 2/"« l ] = 2A/"g [1 - 2 p q 22 (r c )} ( 8 ) 

Here /"q is given by the spontaneous expression Eq.EJ 
and P 22 i. T c ) is the excited state population obtained from 
the time dependent solution of the optical Bloch equa¬ 
tions EL by integration over the coherence time r c of 
the incident x-rays. 

In general, the Bloch equations have to be solved nu¬ 
merically for /? 2 2 (r c ). However, if the coherence time is 
much longer than the Auger decay time (h/T = 1.5 fs for 
Co 2 j? 3 / 2 ), the excited state population reaches an equi¬ 
librium value (see Fig. [2]) and the non-linear contribution 
is given by the analytical expression, 

riiq _ 7qT% ^cohA 3 / (87 r 2 c) _ , , 

/NL ~ /0 (Hv-£ 0 ) 2 + (r/2) 2 + i q r q g coh A 3 /(47r 2 c) 

'-V---' 

p\ 2(00) 

where p^i 00 ) is the equilibrium excited state population 
in the limit r c — > 00 . The non-linear contribution /"^ L 
is seen to have the opposite sign of the spontaneous con¬ 
tribution /g /9 . In the limit of high incident intensity we 
simply have p^i 00 ) = 0-5 and ^-fm, = ~fb' q and the 
sample becomes transparent. 

The increase in excited state population p^ijc) for Co 
L 3 excitation as a function of incident intensity and dif¬ 
ferent coherence times of the incident pulse is shown in 
Fig. [21 It was calculated by numerical solution of the op¬ 
tical Bloch equations, assuming a 20 nm thick Co metal 
film and hui = £ 0 , with the spontaneous peak cross sec¬ 
tion cr§ (Fig. |T](b)), corrected for the r c -dependent en ergy 
bandwidth A£ = 2H\/tt In4/T c of the incident pulse [15|. 
With increasing coherence time r c relative to the core 
hole life time, the Rabi oscillations in the excited state 
population are suppressed. The behavior of p^ 2 { T c) f° r 
our chosen r c values as a function of the total intensity 
per coherent pulse, in units of [mJ/cm 2 /T c ] is shown in 
Fig.[2](b). The stimulated threshold is seen to be lowest 
for coherent pulses in the 3-10 fs range. 

Fig.[3]( a ) shows the effective polarization dependent 
absorption cross section for Co given by Eq.[ 8 ] for three 



FIG. 2: (a) Excited state population p 22 {j c ) as a function 
of the linearly polarized incident intensity Ig for different co¬ 
herence times t c of the incident pulses for the L 3 edge of a 
20 nm thick Co metal film. We assumed resonance excitation, 
heu = £ 0 , linearly polarized light and experimental peak cross 
sections as discussed in the text, (b) p 22 ( 7 c) as a function of 
the incident fluence per coherence time of the pulse. 


values of the incident intensity with f'^ L calculated ac¬ 
cording to Eq.[9] In Fig.[3](b) we illustrate the thick¬ 
ness dependence of the absorption contrast (linear po¬ 
larization) obtained from Eq.[7] for several values of the 
incident intensity. At low intensity the sample trans¬ 
mission decreases with increasing sample thickness d 
due to absorption. However, with increasing intensity, 
the transmitted intensity at large d is seen to decrease 
considerably slower due to stimulated forward scatter¬ 
ing. The magnetic XMCD contrast, plotted in Fig.[3](c), 
first increases with thickness up to a maximum around 
d = 1 /(<r a bsPa) = 17 nm, corresponding to one x-ray ab¬ 
sorption length, before it also decreases. 

Fig-El shows the dependence of the transmitted inten¬ 
sity for the stimulated relative to the spontaneous case 
as a function of the incident intensity, calculated for Co 
metal with d = 20 nm and assuming resonant excitation. 
Both the effective absorption cross section and the trans¬ 
mitted intensity reveal a strong dependence on the inci¬ 
dent intensity, with the spin related XMCD contrast (red 
curve) vanishing faster than the charge related XAS con¬ 
trast (black curve). 

The inset reveals a particularly interesting thickness 
dependence of the transmitted XMCD intensity. For 
a thick sample of 100 nm, the remaining small sponta¬ 
neous XMCD contrast of about 1.5%, which according 
to Fig.[3](c) is greatly diminished by absorption, can ac¬ 
tually be increased by nearly a factor of 5 upon stimula¬ 
tion. 

The stimulated onset in Fig. [I] is predicted to be or¬ 
ders of magnitude lower than for the stimulated effects 
observed before 0 for SASE pulses with an average co¬ 
herence time of about 0.5 fs [16(- This is due to our as¬ 
sumption of self-seeded pulses which besides coherence 
times of ~ 10 fs [3] offer a jitter-free photon energy that 
can be resonantly tuned for maximum cross section. In 
addition, the elastic stimulation threshold is lowered by 
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FIG. 3: (a) Change of the effective polarization dependent ab¬ 
sorption cross section 2A (/o ' 9 + 2/^) for the L 3 resonance in 
Co metal for three incident intensity values, assuming align¬ 
ment of the x-ray propagation direction with the film mag¬ 
netization and the long coherence time limit Eq.[9] The low 
intensity cross sections shown as black curves are nearly iden¬ 
tical to the spontaneous ones shown in red in Fig.[T| (b) De¬ 
pendence of the transmission contrast as a function of sample 
thickness and incident intensity for resonant L 3 excitation of a 
Co metal film due to charge absorption with linearly polarized 
light, according to Eq.JT] (c) Same as (b) for the transmitted 
XMCD contrast. 


the coherent enhancement factor t/coh, which for L 3 ex¬ 
citation of a Co film is ~ 500. 

The dependence of the non-linear contribution on the 
incident intensity, given by Eq.[9] may also be expressed 
in terms of the number of incident photons contained 
in a specific volume. If the volume is chosen to be the 
coherence volume V qk per mode qk : then the associated 
number of photons n q k is referred to as the photon de¬ 
generacy parameter. The incident photons that stimu¬ 
late electronic decays, however, need to be present dur¬ 
ing the total atomic clock decay time H/T which defines 
the sample-specific atomic decay volume Vr. The two 


Photons in atomic coherence volume n r 
IQ’ 2 10' 1 1 10 10 2 



FIG. 4: Dependence of the transmitted intensity according to 
Eq.Q and Eq.[5]in the presence versus absence of stimulation 
for a 20 nm Co metal film as a function of incident intensity. 
The black curve represents the linear polarization or charge 
response [/ t 0 rans]stim/[/t°rans]spon and the red curve is the trans¬ 
mitted XMCD difference intensity [i t “ alis -fttanslstim/KAns - 
-Zttanslspon- The top scale is discussed in the text. The in¬ 
set shows the relative transmitted XMCD contrast for film 
thicknesses of 20 nm (red) and 100 nm (dashed orange) as a 
function of incident intensity. 


coherence volumes are given by 


V qk = A 3 


hjj 

A (Huj) ’ 


V r = A 3 


Tvjj 

2tt 2 T 


( 10 ) 


The number of stimulating photons nr in the volume Vr 
is that in the well-known Kramers-Heisenberg stimulated 
correction term 1 + nr, and it can be expressed in terms 
of the incident polarization dependent intensity and field 
amplitude Eq as, 


n 


9 _ 
r — 


1 A 3 

2 t r 2 c T 


T q 



|-Ezp | 2 


( 11 ) 


On the right we have introduced the zero-point (ZP) 
field Ezr> responsible for spontaneous radiative decays. 
For np = 1 the spontaneous and stimulated scatter¬ 
ing intensities become the same, and the incident field 
Eq is equally effective in driving decays as the ZP field 
I-FzpI 2 = 7 r 2 r/(eoA 3 ) corresponding to one virtual photon 
in the volume Vr. This allows us to equate the intensity 
scale on the bottom of Fig. 0] with the number of photons 
nr on top of the figure, and for our case the ZP field has 
the value Ezr — 4.4 x 10 9 V/m. 
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